Useful transformations: from ion-laser interactions to master equations 
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We show set of transformations which allow to obtain analytic solutions in several quantum-optical 
problems. We start with the ion-laser (time dependent) interaction, continue with the problem of 
a slow atom interacting with a quantized field to end with a master equation that describes losses. 
In all cases it is shown that one may find useful transformations that simplify the problems. 



I. INTRODUCTION 

The purpose of this contribution is to show some transformations that simplify Hamiltonians such that they may 
be treated in an analytic form. We will study various problems, starting with the ion-laser interaction, where we have 
shown that there exists a time dependent transformation that linearizes the Hamiltonian with no approximations [l[ . 
We follow with the interaction of a slow atom with a quantized field, the main result of the paper. In this case the 
atom is affected by the mode shape of the field [|J and we show that the interaction may be simplified as we pass 
from a three-body system to a two-body effective interaction, again with no approximations. Finally we analyze the 
case of the master equation (ME) that describes losses for an anharmonic oscillator [3( , one can transform such a ME 
to obtain a simpler equation where all the superoperators commute, allowing its simple integration 

II. ION-LASER INTERACTION 

We consider the Hamiltonian of a single ion (with unity mass trapped in a harmonic potential in interaction with 
laser light in the (optical) rotating wave approximation [5() 

H = i [p 2 + v 2 (t)x 2 ] + hu 21 A 22 + h\(t) (x, t)A 12 + H.c] , (1) 

A a i, are the operators relating the different electronic transitions (two- level) flip operator for the \b) — > \a) transition 
of frequency uj 2 i, respectively. v(t) is the trap (time dependent) frequency, A the electronic coupling matrix element, 
and E^~\x,t) the negative part of the classical electric field of the driving field. The operators x and p are the 
position and momentum of the centre of mass of the ion. We assume the ion driven by a laser field E^~\x, t) 

E l ~-\x,t) = E Q e- i{ki: - ut \ (2) 

We want to solve the Schrodinger equation 

^™ =%(*)>, (3) 

in order to do this, we make the transformation \(jy) ~ T(t)\£), with [|| 

f (t) = /-^%^ {^) e -i^k &2 (4) 
with pit) a function that obeys the Ermakov equation 

P + v 2 {t)p=\. (5) 

such that we obtain the equation for \4>) 

m^ = Hm)), (6) 
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with the transformed Hamiltonian given by 

n = — ^ if + vl* 2 ) + hw 21 A 22 + nn^le-^p^^-^Au + H.c], (7) 

with fl = XEq. We consider that u> 2 i = u> + 5 where 5 is the so-called detuning. We transform to a frame rotating at 
ui by means of the transformation T w = e ~ lu)tA ^ to obtain the Hamiltonian TL^ = TJHT^ (|</>) — + \4>u>)) 

Hu = M>(t) (n+ M + hSA 22 + hSl(t)[e- l ^+~ a '^A 12 + H.c.}. (8) 

where h = cva with 

a = </ — x + i—= , a' — \ —x ~ i—. (9) 

the annihilation and creation operators respectively. uj(t) = 1/ p 2 is the characteristic frequency of the time dependent 

harmonic oscillator. The time dependent Lamb-Dicke parameter is written as ij(t) — r]Qp(t)^/T^ with t]q = k^J 

where k is the magnitude of the wave vector of the driving field. 

We will consider now the resonant interaction (S = 0). Passing to a frame rotating at the frequency oj(t) we may 
get rid off the harmonic oscillator term in ([5]) to end up with the (time dependent) interaction Hamiltonian 

H(t) = hmie- t{ ^ f * {t ' )dt ' +afelIQ(t ' )dt ' )ri{t) A 12 + H.c.}. (10) 



A. Linearizing the system 

Finally we make the transformation \ip) = R(t)\(f>) with (see Q for the time independent case) 

R(t) = e f(^2l-il2) e -i^(a+a t )(i 22 - J 4 11 ) 

to obtain 



dt 



= h jw(i)n + n(A 22 - i n ) + f | + i[S0(t) - a f /3*(t)]) (ii2 + i 2 i)| |^), (12) 



with /?(*) = _ ir)(t)/2 and we have disregarded the term uj/2 as it would add an overall phase. A method to 
solve J C- like interactions with time dependent parameters has been published by Shen et al. Q . 



B. Many ions 

We generalize the transformation that allows to linearize the ion-laser Hamiltonian 0] in an exact form for two 
different interactions, namely, for the case of many ions in interaction with laser fields [9], and for an ion vibrating in 
two dimensions interacting with a laser field. This linearization has been shown to be important for instance in the 
implementation of fast gates in ion-laser interactions [To| . This is possible, as more regimes may be studied with such 
a linearization: high intensity, low intensity and middle intensity regimes. 

In particular we will show that in the case of many ions, the transformation produces a term which correspond to 
a dipole-dipole interaction. 

Ions in a linear trap interacting with a laser field may be described by the Hamiltonian Q 

H M = vata + - °*3 + Yl %(<7+je i7?i(at+a) + H.c.) (13) 

3 3 

where v is the frequency of the vibration, and a are the creation and annihilation operators of the quantized 
oscillator, 5 is the detuning between the transition frequency of the internal states of the ion, tv eg , and the laser 
frequency ujl, and fl; is the resonant Rabi frequency of the i'th ion in the laser field. The exponentials account for 
the position dependence laser-field and the recoil of the ions upon absorption of a photon. The positions of the ions Xi 
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are replaced by ladder operators kxi = rji(a' + a), where the Lamb-Dicke parameter r\i represents the ration between 
the ionic excursions within the vibrational ground state wavefunction and the wavelenght of the exciting radiation. 
We can linearize the Hamiltonian (|13| via the transformation [3] 



3 



which gives the Hamiltonian 

H M = THT^ 

= va a - 2 2^ ax i + 2^ "j^J + 1 2^ 2 axj + ^ —4- a xj a xk ( 15 ) 

3 3 3 ' j,k 

Here a X j = a+j +<7-j. Note that the transformation, besides linearizing the Hamiltonian, produces an ion-ion (dipole) 
interaction. 

C. Two-dimensional vibration 

An ion vibrating in two-dimensions has the Hamiltonian 

H 2d = v x a\a x + v y a\a y + ia„ + n(a + e"^ +a ^ e ^K+ a *) + H.c.) (16) 

with the transformation 

we can cast the Hamiltonian (11611 into the linearized Hamiltonian 
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T-t2d = v x a\a x + Vya y ay - -a x + fla x + —a x [r\ x {a x - a\) + r\ y (a y - a y )\ (18) 

where we have disregarded a constant term. The Hamiltonian above looks like a two-mode quantized-field-atom 
interaction. 

III. SLOW ATOM INTERACTING WITH A QUANTIZED FIELD 

Here we show how a three body problem may be reduced to a two body problem via a transformation, we treat 
the problem of a slow atom interacting with a quantized field. Because the slowness of the atom, the field mode- 
shape affects the interaction. We can write down the Hamiltonian describing a single two-level atom passing an 
electromagnetic field confined to a cavity. Inn addition to the Jaynes-Cummings Hamiltonian, we have to add the 
energy of the free atom and the spatial variation it feels from the cavity, the Hamiltonian reads 

2 

H = — +uh+ —a z +g(x){aa + + a t cr_), (19) 
on resonance, we can pass to the interaction picture Hamiltonian 

2 

flj= Y+0(aO(aa-+ + aV_), (2°) 

we use the 2x2 notation for the Pauli spin matrices and write the interaction Hamiltonian as (see [J]) 

O'-'+^iJ+T^)* (21) 



where a non unitary transformation T has been used. We define T as 

1 
V 



( ' !' ) <'->*• 
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Note that fft = 1 but f^f = 1 - p g>v with 

|0)<0| 



Ps.v - ( n ln\/nl ) ( 23 ) 



We can use the definitions above to rewrite (jT9|) as 

„2 



(r t T + P 9 ,„)^(T t T + p^)+. 9 (a ; )ft^ ^ / Y T jf (24) 

By noting that Tp gv = we rewrite the above equation as 



where we have used that p 2 g v — p 9l v Finally we factorize the transformation operators in the Hamiltonian above to 
obtain 

2 2 

Hi = f^y + g(x)a x Vn + l)f + V —p g , v (26) 

-v, 2 a 2 

Note that [Tt(^- + g(x)a x y/n + 1)T, ^-p 9 ,t,] = so that the evolution operator for the Hamiltonian above is given by 

Ui(t) = e - if H^+a(x)<y x V : ii+i)ft e ~i^p g ^t ^27) 
to obtain the first exponential we can do Taylor series, and we note that the powers of the argument are simply 



such that 



[^(y + ff(x)a x v^+l)f] & = ^(y + ff(x)<7, V^TT) fc f , fc > 1 (28) 



-if^{S T +g{x)a x s/fi+l)ft _ fj e -i(^-+g(x)<r m VA+l)tf + p gv . (29) 



Note that the evolution operator in (|29p is effectively the interaction of two systems, as it is written in a form in 
which the filed operators commute, unlike the Hamiltonian (119|) . where all the operator involved do not commute. 



IV. MASTER EQUATIONS 

Now we turn our attention to the superoperator solution of master equations for more quantum optical systems, 
namely a dissipative cavity filed with a Kerr medium 0] , master equation describing phase sensitive processes hl| and 
parametric down conversion [12J . Usually these equations are solved by transforming them to Fokker-Planck equations 
hj which are partial differential equations for quasiprobability distribution functions typically the Glauber-Sudarshan 
P-function and the Husimi Q-function. Another usual approach is to solve system-environment problems is through 
the use of Langevin equations, this is stochastic differential equations that are equivalent to the Fokker-Planck equation 

These approaches to the problem makes it usually difficult to apply the solutions to an arbitrary initial field in 
contrast with the superoperator techniques where it is direct the application to an initial wave function. We have 
used this feature in the former Section where we have exploited this fact to obtain reconstruction mechanisms that 
allowed us to obtain information on the state of the quantized electromagnetic field via quasiprobability distribution 
functions. 

A. Kerr medium 



Before applying superoperator methods in the solution of the above equation, let us show how it may be casted 
into a Fokker-Planck equation. In order to do this one writes the density matrix in terms of the Glauber-Sudarshan 
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P-function, p = 1/tt / P(a)\a)(a\d 2 a. Noting that the creation and annihilation operators have the following relations 
with the coherent state density matrix [151 ] 

tf\a)(a\ = ^+aA\a)(a\, (30) 



\a)(a\a 

we can obtain the following correspondence 

dp — > oP(a), 

and 

a*P(a), 



pa^ 



<9 



pa 



a \a)(a\, 



d_ 

da 



P(a), 



(31) 



(32) 



(33) 



In this form, whenever a creation or annihilation operator occurs in the master equation, we can translate this into a 
corresponding operation on the Glauber-Sudarshan P-function. The equation that results is a Fokker-Planck equation 

m 



dP(a,t) 

m 



d 



d 



7 — ol + — — a* + 2"fn 



da 



d 2 



dada* 



P(a,t). 



This equation is equivalent to the stochastic differential equation [14[ 

7a + y/2jfi£(t), 



da 
~dt 



(34) 



(35) 



and the corresponding complex conjugate equation. The quantity is a white noise fluctuating force with the 
following correlation properties 



m = 0, mat)) = s(t - 1 1 ) mat)) = rwr w> = o 



(36) 



Equation (|35[) is usually called a Langevin equation, and also as the " Stratonovich form" of the Fokker-Plank equation 
(see for instance 0, The Fokker-Planck equation can be also obtained from the Kramers- Moyal [H [H 

expansion, a Langevin equation that does not stop after second order derivatives (see for instance [13]). The master 
equation for a Kerr medium in the Markov approximation and interaction picture has the form Q 



dp 
~dt 



-i\[n 2 , p] + 2^dpa) — 70^ dp — ^pd) d. 



(37) 



Milburn and Holmes @ solved this equation by changing it to a partial differential equation for the Q-function and 
for an initial coherent state. We can have a different approach to the solution by again using superoperators. If we 
define 



Yp = -ix[n 2 ,p] 



we rewrite (13"7T) as 



dp 



= (Y + J + L)p, 



where the superoperators J and L are defined as 

Jp = 2^/dpd' , Lp = —^/d)dp — r ypa ! d. 

Now we use the transformation 



(38) 



(39) 



(40) 



p = exp[(f + L)t]p 



(41) 
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to obtain 



with 



^ = exp[-ixRt - 27*] J/3, 



= 2(^,5 — pn). 



(42) 



(43) 



In arriving to equation f4"2"| we have used the formula e yA Be yA — B+y[A, B] +y 2 /2\[A, [A,B]] + . . . for y a parameter 
and A and i? operators. We have also used the commutation relation 

[Y,J]p = 2i x RJp- (44) 

Now it is easy to show that R and J commute, so that we can finally find the solution to equation (|37|) as 

p(t) = e Yt e Lt exp[e J]p(0). (45) 

Note that the above solution may be applied easily to any initial density matrix: 



P(t) 



k,n,m= 



-»yt(» 2 -m 2 )- 7 t(»+m) ■ / 0j± + fc ) ! 



7 1 m I 



-, fc 



2%\{n — m) + 27 



(27)' 
fc! 



\n)(n\ 



(46) 



where /5„ jm (0) are the (Fock) matrix elements of the initial density matrix. 



V. CONCLUSIONS 



We have shown the importance of looking for transformations that simplify Hamiltonians before trying to solve a 
given system. In particular, we have shown how a 3-body system can be reduced to a 2-body system, in the case of a 
slow atom interacting with a quantized field. We have given the most complete solutions for the ion-laser interaction 
in several cases: time dependent case, several ions, ions vibrating in two dimensions. 
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